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ABSTRACT. We analyze the optimal dividend payment problem in the dual model under constant trans- 
action costs. We show, for a general spectrally positive Levy process, an optimal strategy is given by 
a (ex, C2)-policy that brings the surplus process down to c\ whenever it reaches or exceeds C2 for some 
< ci < C2. The value function is succinctly expressed in terms of the scale function. A series of nu- 
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merical examples are provided to confirm the analytical results and to demonstrate the convergence to the 
no-transaction cost case, which was recently solved by Bayraktar et al. [8]. 
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1. Introduction 

We solve the optimal dividend problem under fixed transaction costs in the so-called dual model, in 
which the surplus of a company is driven by a Levy process with positive jumps (spectrally positive 
i/-) ; Levy process). This is an appropriate model for a company driven by inventions or discoveries. The case 
without transaction costs has recently been well-studied; see [3], [7], [2], and [4]. In particular, in [8], we 
show the optimality of a barrier strategy (reflected Levy process) for a general spectrally positive Levy 
process of bounded or unbounded variation. 

A strategy is assumed to be in the form of impulse control; whenever dividends are accrued, a constant 
transaction cost (3 > is incurred. As opposed to the barrier strategy that is typically optimal for the 
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no-transaction cost case, we shall pursue the optimality of the so-called (ci, c 2 )-policy that brings the 
surplus process down to c\ whenever it reaches or exceeds c 2 for some < c\ < c 2 < oo. While, as 
in [15], an optimal strategy may not lie in the set of (ci, c 2 )-policies for the spectrally negative Levy 
case, we shall show that it is indeed so in the dual model for any choice of underlying spectrally positive 
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Levy process. As a related work, we refer the reader to a compound Poisson dual model by [20] where 
transaction costs are incurred for capital injections. 

Following [8], we take advantage of the fluctuation theory for the spectrally positive Levy process 
(see e.g. [9], [11] and [14]). The expected net present value (NPV) of dividends (minus transaction 
costs) under a (c 1; c 2 ) -policy until ruin is first written in terms of the scale function. We then show the 
existence of the maximizers < c* < c 2 < oo that satisfy the continuous fit (resp. smooth fit) at c 2 when 
the surplus process is of bounded (resp. unbounded) variation and that the derivative at c* is one when 
c* > and is less than or equal to one when c\ = 0. These properties are used to verify the optimality of 
the (c*, c 2 )-policy. 

In order to evaluate the analytical results and to examine the connection with the no-transaction cost 
case developed by [8], we conduct a series of numerical experiments using Levy processes with positive 
i.i.d. phase-type jumps with or without Brownian motion [1]. We shall confirm the existence of the 
maximizers < c\ < c 2 < oo and examine the shape of the value function at c\ and c 2 . We further 
compute for a sequence of unit transaction costs and confirm that, as /3 4- 0, the value function as well as 
c\ and c 2 converge to the ones obtained for the no-transaction cost case in [8]. 

The rest of the paper is organized as follows. Section 2 gives a mathematical model of the problem. 
In Section 3, we compute the expected NPV of dividends under the (ci, c 2 )-policy via the scale function. 
Section 4 shows the existence of < c\ < c 2 < oo that maximize the expected NPV over c\ and c 2 . 
Section 5 verifies the optimality of the (c*, c* 2 ) -policy. We conclude the paper with numerical results in 
Section 6. 

2. MATHEMATICAL FORMULATION 

We will denote the surplus of a company by a spectrally positive Levy process X = {X t ; t > 0} 
whose Laplace exponent is given by 

(2.1) ifj(s) :=logE [e~ sXl ] =cs + \a 2 s 2 + [ {e~ sz - 1 + szl {0<z<1} )u{dz), s eR 

* J{0,oo) 

where v is a Levy measure with the support (0, oo) that satisfies the integrability condition J^ Q ^(1 A 
z 2 )u(dz) < oo. It has paths of bounded variation if and only if a = and f, Q ^ z v(dz) < oo; in this 
case, we write (2.1) as 

ip{s)=T)s+ f (e~ sz - l)v(dz), seR 

i(0,oo) 

with := c + f. Ql . zu(dz). We exclude the trivial case in which X is a subordinator (i.e., X has 
monotone paths a.s.). This assumption implies that o > when X is of bounded variation. 

Let F x be the conditional probability under which X = x (also let P = P )> an d let F := {J^ : t > 0} 
be the filtration generated by X. Using this, the drift of X is given by 

(2.2) /i := E[Xx] = -^'(0+). 
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In order to make sure the problem is non-trivial and well-defined, we assume throughout the paper that 
this is finite. 

Assumption 2.1. We assume that ji G (—00, 00). 

A (dividend) strategy ix := {LJ; t > 0} is given by a nondecreasing, right-continuous and F-adapted 
pure jump process starting at zero in the form I7 t = J2o< s <t ^1 w ^ tn ^L t = L t — L t -, t > 0. 
Corresponding to every strategy n, we associate a controlled surplus process U n = {U£ : t > 0}, which 
is defined by 

U? := X t -LI, t> 0. 

The time of ruin is defined to be 



a 77 := inf {t > : < 0} . 
A lump-sum payment cannot be more than the available funds and hence it is required that 



(2.3) 



AL* < Ut_, t < a* a.s. 



Let IT be the set of all admissible strategies satisfying (2.3). The problem is to compute, for q > 0, the 
expected NPV of dividends until ruin 



VAX) 



^° 0<s<t 



x > 0, 



where (3 > is the unit transaction cost, and to obtain an admissible strategy that maximizes it, if such a 
strategy exists. Hence the (optimal) value function is written as 



(2.4) 



v\x) 



sup V n (X ) 
Tren 



x > 0. 



3. The (ci, c 2 ) -policy 



We aim to prove that a (c\, c^-policy is optimal for some c* 2 > c\ > 0. For c 2 > C\ > 0, a (ci, c 2 )- 
policy, 7r Cl|C2 := {-^j 1 ^ 2 ; t > 0}, brings the level of the controlled surplus process U Cl ' C2 := X — L Cl ' C2 
down to ci whenever it reaches or exceeds c 2 . Let us define the corresponding expected NPV of dividends 



as 



(3.1) 



J Cl,C 2 



» := E a 



0<s<t 



{AL?' C2 >0} 



x > 0, 



where cr CljC2 := inf {t > : f/ t Cl ' 62 < 0} is the corresponding ruin time. In this section, we shall express 
these in terms of the scale function. 
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3.1. Scale functions. Fix q > 0. For any spectrally positive Levy process, there exists a function called 
the q-scale function 

W {q) : R i — y [0,oo), 

which is zero on (— oo,0), continuous and strictly increasing on [0, oo), and is characterized by the 
Laplace transform: 

1 



"W {q \x)dx 



ip(s) -q 



s > 



where 



$(g) := sup{A > : ip(X) = q). 

Here, the Laplace exponent ip in (2.1) is known to be zero at the origin and convex on R + ; therefore <3>(g) 
is well-defined and is strictly positive as q > 0. We also define, for i6R, 



W (q \x) :-- 



Z {q) \ 



x 



W^(y)dy, 



1 + qW (q \x) 



Z {9 \x):= Z^(z)dz = x + q / W iq) (w)dwdz. 
Jo Jo Jo 

Notice that because is uniformly zero on the negative half line, we have 

(3.2) Z {q) (x) = 1 and Z {q) (x) = x, x < 0. 
Let us define the first down- and up-crossing times, respectively, by 

(3.3) r~ := inf {t > : X t < a} and r 6 + := inf {t > : X t > b} , a,b e 
Then we have for any b > 

W^ib-x) 



(3.4) 



E, 



E, 



e -ir„ I 



e- qT * 1 



W>ro} 



W<ro} 



WM(b) 



Z^ q \b 



x) 



Z^ q \b) 



W( q \b-x) 



WM(b) ' 

Fix a > and define ip a {-) as the Laplace exponent of X under P a with the change of measure 



(3.5) 



dP 



exp(aX t — ip(a)t), t > 0; 



as in page 213 of [14], for all s > —a, 

• i 



aa + c 



me 



o 



l)u(du))s + -a 2 s 2 + \ 



1 + snl {0<u<1} )e au u(du). 
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If W^ q) and Z { a q) are the scale functions associated with X under P a (or equivalently with i/) a (-)). Then, 
by Lemma 8.4 of [14], 



(3.6) 



W (Q-^(a)) ( x ) = e -ax w ( q ) xe 



which is well-defined even for q < ift(a) by Lemmas 8.3 and 8.5 of [14]. 

Remark 3.1. (1) If X is of unbounded variation, it is known that W^ q ' is C 1 (0, oo); see, e.g., Chan 
et al. [10]. Hence, 

(a) Z^ q > is C 1 (0, oo) and C° (R) for the bounded variation case, while it is C 2 (0, oo) and C 1 { 
for the unbounded variation case, and 

(b) isC 2 (0, oo) andC 1 (M.) for the bounded variation case, while it is C 3 (0, oo) andC 2 ' 
for the unbounded variation case. 

(2) Regarding the asymptotic behavior near zero, we have that 



(3.7) 



W(«)(0) 



0, if X is of unbounded variation, 
|, if X is of bounded variation, 



and 



(3.8) 



W {q) '(0+) := \imW {qY (x) 

x\.0 



OO, 

q+^(0, oo) 



(3) As in (8.18) and Lemma 8.2 o/[14], 

W®\y) W^'(x) 



< 



W(i){y) ~ WW 



x 



if a > 0, 

if a = and z/(0, oo) = oo, 
if X is compound Poisson. 



y > x > 0, 



where W^' is understood as the right-derivative if it is not differentiable. In all cases, (x—)> 
W^' (x+) for all x > 0. 

3.2. The expected NPV of dividends for the (ci, c 2 )-policy. Now we obtain (3.1) using the scale func- 
tion. By the strong Markov property, it must satisfy, for every < x < c 2 and < c\ < c 2 , 



(3.9) 



v c u c 2 (x) = ^x 



-QTc 



= 1 



y ci,c 2 1 



where v cltC2 := v CuC2 (ci). Solving for x = c\, we have 



(3.10) 



-ei,C2 



E C1 


e-^l {r _ 




-0)\ 


1-E C1 


e "" Tc+2l K + 2 <nr}. 





< ci < c 2 . 



We will express these quantities in terms of the scale functions using (3.4) and the lemma below. 
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Lemma 3.1. For < x < c, 



E, 



-K+<r ( n 



-R {q \c-x) + ( c- ^ ) 



C JW(?)( C )_#( ( 

9/ 



W^)( c ) ' 



9 



Proof. Suppose X = —X and let be the law under which X = x. The stopping times r are defined 
analogously to (3.3). Because X T - = a.s. on {t ~ < oo} under F x for x > 0, 



E, 



E,, 



e -9(r+Ar ) X + 



t c Ar„ 



-E_ 



-?(t- +At - c )X~ 



To compute this last expression, we consider the moment generating function: 



F(x; c,v) := E_ x 



By changing measure using (3.5), and then applying (3.4) and (3.6), we obtain 



F(x;c,v) = e- vx E 



g — (c — x) 



e~ vx E v 



W^ {v) \c-x) 
W^ (v)) (c) 



+ Z^-^\c - x) - Z^~^\c) 



e -(q-iP(v))(T+ at_ (c _ x) ) 

W^ {v) \c-x) 



Z^)\ C -x)-(q- ^W^^ic] Wv 



(q-i[)(v)) 



-vxrxlil-^)) ( . _n TrM-^)) , _ % (C - X) 



e~ vx w: 



[c-x)-w: 



WM(c) 



Since we assumed that fi = —ip , (0+) G (— oo, oo), we can write 



E, 



-gr c 



-{t+<t-} X t+ 



hm — F (rr; c, w) 



(3.11) 



x — fl 



W (c — x) — W cr 



— q lim 



'c- X )-Wf {v)) (c) W{q){c - x) 
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Using e- vx W { r^ V)) (c - x) = J^~ x e- viv+x) W iq \y)dy (which follows from (3.6)) and integration by 
parts, we obtain 



hm — e VV„ c — x = — xW (c — x 

v^o dv 



yWM(y)dy 



o 

PC — X PC~~X 

-xW iq \c-x)-(c-x)W iq \c-x)+ / W iq \y)dy = -cW i9 \c-x)+ / W {q \y)dy 

Jo Jo 



X 



cZ( 9 \c-x)+Z l9 \c-x) 



Again using integration by parts, it follows that \im. v -+o(dW^ ^ ^ (c)/dv) equals 



- lim / ye- vy W {q \y)dy = - [ yW {q \y)dy = -cW (q) (c) + [ W (q) (y)dy = - \-cZ {q \c) + Z {q \c) . 
v ^°Jo ' Jo Jo Q L J 



Substituting these in (3.1 1), we get that E 2 



e i K+<r -} A r^ 



equals 



x — 



Z^{c-x)-l-{Z (q \c)-l) 



W {q \c-x) 



-x + cZ {q \c -x)- Z {q \c -x)- (cZ®(c) - Z w (c 



W^jc-x) 



.a 

Q 



Z {q \c-x) - Z {q \c) 



+ cZ^ q \c -x)- R {q \c -x)+ (-cZ®(c) + R {q) (c)) 



W^(c-x) 
f(«)(c) 



-R {q \c-x) + c- - Z^ q \c — x) 



c-^j Z^\c)-R {q \c) 



x) 



WM(c) ' 



as desired. 



□ 



By this lemma, (3.4) and (3.10), we can write 



(3.12) 



J Ci,C2 



f(ci,c 2 ) 

g(c u c 2 ) 



< ci < c 2 , 
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where 



/(ci,c 2 ) 



-R {q \c 2 - Cl )+ (<*-£) Z^(c 2 - Cl ) 



c 2 -^j Z^\c 2 )-R^\c 2 ) 



W^(c 2 - Cl ] 



(3.13) 



(ci + /3) 



c 2 — C\ 



Z^\c 2 ) 



WM(c 2 ) 



-^(c 2 - Cl )+^ 2 - Cl -/3-^ Z<«> 



W^(c 2 ) 

c 2 - CI', 



c 2 - Ci 



0-^ Z«(cs)-i2 (ff) (ca) 
9/ 



^ (g) (c 2 - Cl ) 
^(«>(c 2 ) 



and 
(3.14) 



9{ci,c 2 ) 



c 2 - Ci 



W(9)(C2) 



4. Candidate strategies 



Using the results in the previous section, we now have an analytical expression for (3 . 1 ) or equivalently 
(3.9). For < x < c 2 and < c\ < c 2 , this expression reduces to 



v ci ,c 2 (x) = -R [q \c 2 -x) + 



c 2 



(4.1) 



c 2 - 



Z^\c 2 )-R^\c 2 ) 



+ ( V d,c 2 - Ci - f$) 



z {q \ 



c 2 - X) 



c 2 - X) 

W^i^-x) 

W^{c 2 ) 
z{q){ ^W^(c 2 -x) 



c 2 J 



where 
(4.2) 

For a; > c 2 , we have 
(4.3) 



-R^(c 2 -x) + 7 (ci, c 2 )Z^(c 2 - x) - G{a, c 2 



7(ci, c 2 ) := f Cl . C2 + c 2 - ci - /3 - -, 

q 

G(c u c 2 ) :=j( Cl ,c 2 )zM(c 2 )-RM(c 2 ) 



v CuC2 (x) = x - Ci - (3 + v c 



W(i)(c 2 ) 

W^(c 2 -x) 



W^(c 2 ) 



In view of (4.3), a necessary condition for a (c±, c 2 )-policy to be optimal is that c\ and c 2 maximize 
c\. In this section, we first obtain the first-order conditions by computing its partial derivatives 



J Cl,C2 



with respect to c\ and c 2 and then show the existence of finite- valued maximizers. In the rest of the paper, 



OPTIMAL DIVIDENDS IN THE DUAL MODEL UNDER TRANSACTION COSTS 



the derivative is understood as the right-derivative when the scale function fails to be differentiable 

on (0, oo). 







4.1. First-order conditions. 
Lemma 4.1. For every < c\ < c 2 , 

-K- {Vci,c 2 ~ Ci) = ^~ 

oc 2 dc 2 
Proof. Differentiating (3.13), we obtain 



J Cl,C 2 



G{c u c 2 ) d W®(C2- C1 ] 

' g(c u c 2 ) dc 2 W^(c 2 ) 



d 



dc' 



/(ci,c 2 ) 



+ Z^(c 2 - Cl ) + (c 2 - Cl - p - ^) qW^(c 2 - Cl ) 



c 2 ~ Cl 



P-£) <lW {q \c 2 ) + Z^\c 2 ) - Z^\c 2 ) 



W^(c 2 - Cl ] 



Z^(C 2 )-RM(C 2 ) 

qJ 



c 2 - Cl -(3-^)Z^(c 2 )-R^(c 2 
q 



W(i)(c 2 ) 
d W^(c 2 - Cl ) 
_ 8c~ 2 W^(c 2 ) 

d W {q) (c 2 - Cl ) 



On the other hand, differentiating (3.14) yields 



dc 2 WM(c 2 ) ' 

d W^(c 2 - Cl ) 



c 2 



d W^{c 2 -d) 



dc 2 WM(c 2 ) 
Using the last two equations along with (3.12), we have 

d d d 

g(ci,c 2 )—v CU2 = ^-/(ci,c 2 ) - v CUC2 —g(c u c 2 ) 



c 2 - Ci 



P-£) Z^\c 2 )-R^ 



-G(ci,c 2 ) 



dc~ 2 W^{c 2 ) 



c 2 ; 



d W^(c 2 - Cl ) 
_ dc~ 2 W^{c 2 ) 



J Cl,C 2 



C 2 



9 w( q \c 2 - Cl ] 
dc~ 2 W("\c 2 ) 



□ 



Lemma 4.2. For < c\ < c 2 , 



d Cl [Vci ' C2 

where 



d f f(c h c 2 ) - dg(c u c 2 ] 



dci 



9{Ci,c 2 ) 



9{ci,c 2 ) L 



H(ci,c 2 ) + G(ct,c 2 



H( Cl ,c 2 ) := q j( Cl ,c 2 )W^(c 2 - c x ) - W (q \c 2 - cy 



c 2 - Ci 



W(i\c 2 ) 
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Proof. By (3.13) and (3.14), 

f( Cl ,c 2 ) - c ig ( Cl ,c 2 ) = -R^{c 2 - d) + (c 2 - Cl -P - Z^(c 2 - ci) 



c 1 ^(c 2 -c 1 )-c 1 Z«(c 2 )55^ 



and 



hence its derivative equal 

9ci 



-ci + 

= -i? (9) (c 2 - d) - Cl +(c 2 -p-^) Z^(c 2 - Cl ) 



[/(ci, c 2 ) - d#(d, c 2 )] = gW (?) (c2 - ex 



- c 2 



+ 



i x d ( f(c u c 2 



C 2 K l 

= qW^(c 2 - Cl ) - ZW(c /^) Cl) ^d by (3.12), 

5 



Because dg(ci, c 2 )/dci 

-d#(d,c 2 )\ 9 . , . o» 

= ^-[/(Cl,C 2 ) -Ciflf(Ci,C2)J - (t' Cl , C2 ~ C 1 )—g{C h C 2 



(c 2 - /3 - ^ gW^(c 2 - d) + (c 2 - /3 - ^ zW(ca) - it^c. 



= giy (9) (c 2 - Cl 

- (e M - ci) [gW^(c 2 - Cl ) - ^ (9 Hc 2 ) ^ ( ^ ( ( g ) 2 ( ~ Cl 

W(9)'f c „ _ Cl ) 
= -F( Cl)C2 ) + G ( Cl , C2 )_^^. 



V^'(c 2 -d) 
V^)(c 2 ) 



□ 



Remark 4.1. The first-order conditions obtained above are for (4.3). However, these are in fact the same 
for (4.1) for any < x < c 2 . Differentiating the first equality of (4.1), 



oci ac\ 



c 2 - a;) 



C2 



W^(c 2 -x) 
W^{c 2 ) 



< ci < c 2 , 



whose sign is the same as that ofd(v CltC2 — d)/ dc\ thanks to (3.4) which guarantees that the expression 
inside the square brackets is positive. Moreover, by differentiating (4.1) and by Lemma 4.1, for < c\ < 
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C 2 , 







do C1,C2 



-C7(ci,c 2 ) 
-G(c u c 2 ) I 



x) 



dc 2 WM(c 2 ) 

Z®(cn-x) 



X) 



W^(c 2 -x) 



fi , (ci,c 2 ) 



W(«)(c 2 ) 

d W^(c 2 -x 
dc~ 2 W^{c 2 ) 



_d_ 

dc-2 



whose sign is the same as that ofdv CltC2 /dc 2 due to item (3) of Remark 3.1. 

4.2. Existence and some properties of maximizers. Now we are ready to show that the maximizers of 
V C1 ,C2 ~ c i exist. We will also describe equations that can be used to identify these points. 

Lemma 4.3. We have sup 0<Cl<C2 (v CliC2 — ci) = sup 0<Cl<C2<c (w Cl)C2 — ci) for sufficiently large C < oo. 

Proof. By Lemma 3.1, for any c 2 > C\ > 0, 



E 



e _97 M, 



K+2<T - } (X Tc+2 - f>)\ = -RM(c 2 - Cl ) + ( c 2 - ft - t \ Z ^\c 2 - Cl ) 



Q 



c 2 -(3-^)z^(c 2 ) 



W^(c 2 - Cl ] 
W^(c 2 ) 



and hence 
(4.5) 



_d_ 

dci 



E, 



c 2 -P 



)Z^(c 2 )-R^(c 2 ) 



dc~ 2 W^{c 2 ) 



A(c)W®(c c ) ( W(qY{C2 ~ Cl) WiqY(C2) 



where A(c) := (c-/3-^j - o0.lt follows from Exercise 8.5 of [14] and Proposition 

2 of [5] that Zfo)(c)/WW(c) g/$(g) e (0, oo) and #(c)/f(«'(c) -»■ g/$(g) 2 e (0, oo) as c t oo, 
respectively. As a result, A(c) t °° an d hence there exists B < oo such that 



(4.6) 



A{c) > 0, c > B. 



Now because ^gfg-W^ > by Remark 3.1(3), we have <9E Cl [ e -*4 l {r + <T - } (X T + - /?)] /<9c 2 < 
for any c 2 > C\ > B. Hence for any fixed c\ > B, 



sup E Cl 

C2X2>Ci 



lim E, 



Cl 
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Now by the definition of v CljC2 as in (3.10), for any fixed c\ > B, 



sup (v Cl ,c 2 - ci) = sup 

C2:C2>ci C2-C2>c\ 



-ci + E C1 




l-E Cl 


6 2i {Tc+<T -}_ 





< sup 

c 2 :c 2 >ci 



-ci + sup C2:C2>Cl E C1 


e " 9Tc+2l {4<^}(^ 2 -^ 


1-E C1 


p~<l T Cn 1 , 

e 2i {r+<r ( 7}_ 





< sup 

c 2 :c 2 >ci X - E Cl 



-P- A(ci)W®(0) 



e -qr C2 I 



K + 2 < r 1 



which is negative by (4.6). On the other hand, because c\ — and c 2 > attain u CliC2 — c\ — 0, we have 
(4.7) 



sup (u cl]C2 - c iy 

(ci,C2):C2>ci>0 



sup (_U C1)C2 - ClJ 

(ci,c 2 ):C2>ci>0,ci<S 



Now fix ci < B and c 2 > £> + 5 for any 5 > 0. Then 



Cl 



-ci + E C1 






-B-/3 


1-E C1 


6 2i K + 2 <-0"}. 


" 1-Eb 


t T S+<5 <7 J_ 



=: M > -oo. 



By Lemma 4.1 and Remark 3.1(3), 



d^ VcuC2 



G( Cl ,c 2 ) d wMfa-d] 
' g(c u c 2 ) dc 2 WM(c 2 ) 



< 



~W,c 2 + C 2 ~ Ci - P 



Z^(c 2 ) _ R( q) {c 2 ) \ W^{c 2 - Cl ) ( WW{C2 - ci) _ W^'(c 2 ) 



W^(c 2 ) W^(c 2 ] 



9{cuc 2/ 



c 2 - Ci 



w/(<?) 



C2, 



M + c 2 - /3 



q 



Z^(c 2 ) R^ q \c 2 ) \ W^\c 2 - Cl ) f W^'(c 2 - Cl ) _ wW(c 2 

~ W(i){c 2 )J g( Cl ,c 2 ) \W^(c 2 - Cl ) ~ WM{c 2 ) ) ' 



Using Remark 3.1(3) and the fact that (M + c 2 - ft 
exists a sufficiently large constant C such that 



-) „. fo A ; , — ... fo A ; >• oo, it follows that there 



sup — u 

ci<B,c 2 >C CC 2 



"Ci,C2 ^ 0. 

OC-) 



Combining the last inequality with (4.7) completes the proof. 
Lemma 4.4. Fix any c\ > 0, lim C2 ^ cl G(ci, c 2 ) < 0. 
Proof. We have 



□ 



/ \ C2\,C\ 

7(ci,c 2 ) > 



— oo, 

Z^(ci)- 1 



iyw (ci) 

' IVW(0) 



/3+i? (9) (ci: 



if X is of unbounded variation, 
if X is of bounded variation. 
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When X is of unbounded variation lim C2 ^ Cl G(ci, c 2 ) = — oo while when X is of bounded variation, by 
(4.2), lim C2 ; Cl C7( Cl , c 2 ) = -|^)W/3 < 0. 

□ 

This lemma, together with Lemma 4. 1 and Remark 3 . 1 (3), implies that, for any fixed c\ > 0, dv Cl jC2 / <9c 2 
is negative near c\\ consequently there exist v CltC1 := lim C2 j, cl v Cl:C2 (which can be shown to be — oo 
when X is of unbounded variation). Because v cliC2 — c\ is continuous and we have a compact domain 
{(ci, c 2 ) : < Ci < c 2 , < c 2 < C} for large C by Lemma 4.3, we have a maximum. Furthermore, 
Lemmas 4.1 and 4.4 show that if c\ and c 2 maximize v Cl ,c 2 — c i> it must hold that c 2 is away from c\. 

Lemma 4.5. Suppose c\ and c 2 maximize v cliC2 —ci. ThenG(ci,c 2 ) = andHic^, c 2 ) > 0. In particular, 
ifci > 0, we must have H(ci,C2) = 0. 

Proof. By Lemmas 4.3 and 4.4, c 2 £ (ci, oo). Hence, by Lemma 4.1, we must have G(ci, c 2 ) = 0. On 
the other hand, by Lemma 4.2, 

9 #(ci,c 2 ) 

^— (^ci,c 2 - Cl) = ? r-. 

dci g{ci, c 2 ) 

If if(ci, c 2 ) < 0, the derivative is positive and it violates the assumed optimality. In particular, if c\ £ 
(0, c 2 ), then the derivative must vanish and hence H{c\, c 2 ) =0. □ 

Combining the above arguments, we arrive at the following proposition. 

Proposition 4.1. There exist (ci, c 2 ) that maximize v cifi2 — c\ and satisfy the following two properties. 

(1) < c 2 < oo and G(ci, c 2 ) = 0; 

(2) either < c\ < c 2 with H(ci, c 2 ) = 0, or c\ = with H(ci, c 2 ) > 0. 

Remark 4.2. Suppose c\ and c 2 are such that H(c\, c 2 ) > and G(c\, c 2 ) = 0. Then, 7(01, c 2 ) > 0. To 
5ee w/zv ?/zz5 z'5 50, />v Lemma 4.4, G(ci, c 2 ) = implies that c\ < c 2 an<i, together with H(ci, c 2 ) > 0, 
we have 7(01, c 2 ) > W W (c 2 - Ci)/W^(c 2 - c t ) > 0. 

5. VERIFICATION OF OPTIMALITY 
By Proposition 4.1, there exist < c\ < c* 2 < 00 such that G(c{, c* 2 ) = and either 

Case 1: c\ > with #(c*, c 2 ) = 0, or 
Case 2: c\ = with if(c*, c£) > 0. 

We will show that such a (c*, c 2 ) -policy describes an optimal policy (and as a result the conditions written 
in terms of H and G are both necessary and sufficient for (c*, c 2 ) to be optimal.) Propositions 5.1 and 
5.2 will play a key role. 
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By substituting G(c{, c* 2 ) = in (4.1), 



V ct,c*{X) 



-R^(c* -x)+ 7(cJ, c* 2 )Z^(c* 2 -x), < x < c*, 
x-c\ -(3 + v chc *, x>c* 2 . 



In fact, by (3.2) and by the definition of 7(0^, c 2 ) as in (4.2), we can write for any x > 0, 

(5.1) r, c ., c .(x) = -RW(c* 2 -x) + 1 (c* 1 ,c* 2 )ZM(c* 2 -x). 
It is clear that it is continuous at c 2 . Regarding its differentiability, we have 

(5.2) v' clA {x) = Z^(c* 2 -x)- 7 (cJ, c* 2 )qW^(c* 2 - x), 
whose limit equals 

(5.3) t^te-) = 1-7(^,4)^(0). 

Because v' c * c *(c 2 +) = 1, the differentiability at c 2 is satisfied if and only if X is of unbounded variation 
by (3.7). We summarize these observations in the lemma below. 

Lemma 5.1 (smoothness at c 2 ). The function v c * c *(-) is continuous (resp. dijferentiable) at c 2 when X 
is of bounded (resp. unbounded) variation. 

Remark 5.1. Differentiating (5 .2) further, 

(5.4) v^(x) = -qW^(c* -x) + jKc* 2 )qWW(c* 2 - x), 
for a.e. x G (0, c 2 ) and its limit as x t c 2 equals 

(5.5) «W<$") = -qW^(0) + 7(c*,c;)^^'(0+). 

These results on the second derivative are used in deriving Propositions 5.1 and 5.2 below. 

By Remark 3.1(1) and Lemma 5.1, the function v c * tC * is C°(0, 00) andC 1 ((0, oo)\{c2}) (resp. C^O, 00) 
and C 2 ((0, oo)\{c2})) when X is of bounded (resp. unbounded) variation. 

Let £ be the infinitesimal generator associated with the process X applied to a sufficiently smooth 
function / 

Cf(x) := -cf'(x) + -a 2 f\x) + J [f(x + z) - f(x) - f'(x)zl {0<z<1} ] u(dz). 

Here £v c * iC *(-) makes sense anywhere on (0, oo)\{c2}. 

Proposition 5.1. (1) (£ — q)v c *, c *(x) = Ofor < x < c 2 , 
(2) (£ - q)v c * A (x) < Oforx > c* 2 . 



OPTIMAL DIVIDENDS IN THE DUAL MODEL UNDER TRANSACTION COSTS 15 

Proof. (1) By Proposition 2 of [5] and as in the proof of Theorem 8.10 of [14], the processes 

c 2 c 2 

are martingales. Thanks to the smoothness of and RS q > on (0, c%) (see Remark 3.1(1)), we obtain 
(C — q)R^ (y) = (£ — q)Z^ (y) = for any < y < c 2 . This step is similar to the proof of Theorem 
2.1 in [8]. This implies claim (1) in view of (5.1). 

(2) Suppose X is of bounded variation. By (5.3) and Remarks 3.1(2) and 4.2, 

Because a = 0, o > and v c * c *(-) is continuous at c?;, (1) implies (£ — g)^ c *, c *(c2+) < 0. Because, on 
(c 2 , oo), £t> c * iC * is a constant and gt> c * jC * is increasing in view of (5.1), claim (2) follows for the bounded 
variation case. 

Suppose X is of unbounded variation. By (5.5) and Remarks 3.1(2) and 4.2, 

v'^(4-) = l(clcl)qWW(0+) > = v'^ c ,(4+). 

Because v c * iC * is differentiable at c 2 , we must have (£ — q)v c * tC *(c 2 +) < 0. Again, because (£ — q)v c * jC * 
is decreasing on (c 2 , oo), (2) is proved for the unbounded variation case as well. □ 

Proposition 5.2. For any x > y > 0, it holds that v c * tC * (x) — v c * c * (y) > x — y — (3. 

In order to show this proposition, we take advantage of the slope of v c * c * at c\. By (5.2), 

< lA {^) = Zi9) ^ - c *) - l(cl,c;)qW^(c* 2 - c\) = 1 - H(4,4). 
When c\ = 0, the derivative is understood as the right-derivative. Hence we arrive at the following. 

Lemma 5.2 (slope at c*). For both Cases 1 and 2, v' c * c * (c*+) < 1. In particular, for Case 1, v' c * c » (c*) = 
1. 

Lemma 5.3. For any x G (0, oo)\{c2}, v' c ,^{x) < 1 if and only if x G (c*,^). 

Proof. Because v' c * c *(x) = 1 on (c 2 , oo), we shall focus on x G (0, c 2 ). Rewriting (5.4), 

(5.6) v': hc ,(x) = -qW^(c* 2 -x)J(x), < x < c* 2 , 

where J(x) :— 1 — 7(0* , c%) ~^ . By Remarks 3.1(3) and 4.2, J(-) is decreasing on (0, c 2 ), and 

hence there exists a unique level c G [0, c 2 ] such that (5.6) is negative if and only if x < c. In other 
words, there are three possibilities 

(i) ^c*, c * is strictly concave on (0, c 2 ) or 

(ii) f c *,c* is strictly concave on (0, c) and strictly convex on (c, C2), 

(iii) f c *,cj is strictly convex on (0, c 2 ). 



16 



E. BAYRAKTAR, A. E. KYPRIANOU, AND K. YAMAZAKI 



Case 1: Suppose c\ > with H{c\, c 2 ) = 0. By Lemma 5.2, (5.3) and Remark 4.2, 
(5-7) <^(cH<l = <*, c *(ct). 

Therefore we can safely rule out (iii) and we must have either (i) or (ii) with c\ < c < c 2 . For (i) (thus 
v' c * iC * is decreasing on (0, c 2 )), given x G (0, c* 2 ), v' c * c *(x) < 1 if and only if x G (c*, c* 2 ). Now suppose 
(ii) with c* < c < c 2 . Then by the concavity on (0, c) and 1 = v' c * )C * (c[), we have v' c * jC * > 1 on (0, cl) and 
v f c ^ c * < 1 on (cl,c). Forx G (c, C2), by the convexity on (c,c 2 ) and (5.7), 1 > v' c , c ,(c 2 -) > v' c * >c *(x). 

Case 2: Suppose c* = with if (0, C2) > 0. In view of (5.2) and the definition of H(c\, c 2 ), we must 
have that v' c * (0+) < 1. This together with v' Q c , (c 2 —) < 1 shows that v' c , (x) < 1 on (0, c^) for any of 
(i), (ii) and (iii). □ 

By Lemma 5.3, 

inf [v c *, c *(x) - v clA (y) - (x - y - P)] = v c * jC *(c* 2 ) - v c * A {c\) - (c* 2 - cl - (3) = 0, 

and as a result the claim in Proposition 5.2 follows immediately. 
Next, we will verify the optimality of the (c*, c 2 ) -policy. 

Theorem 5.1. We have v c * tC *(x) = sup^n v n (x)for every x > and the (c{, c 2 )-policy is optimal. 

Proof. Here we only provide a sketch of a proof since it is similar to that of Lemma 6 of [15]. To verify 
the optimality of (cj, c 2 ) we only need to show that v c * >c *(x) > v n (x), x > 0, for all n G II. But this 
result follows from applying the Ito formula to v c * yC * (U^) for an arbitrary tt G II, using Propositions 5.1 
and 5.2 and then passing to the limit using Fatou's lemma. Here one should be careful in applying the 
Ito formula since the value function v c * tC * may not be smooth enough at c 2 to apply the usual version. 
When X of unbounded variation, we use Theorem 3.2 of [17], which shows that the smooth fit principle 
(which we proved in Lemma 5.1) is enough to kill the local time terms that might accumulate around c 2 ; 
see also Theorem IV.71 of [18], or Exercise 3.6.24 of [13]. On the other hand, when X is of bounded 
variation recall from Lemma 5.1 that the value function is only continuous. However, in this case we do 
not need the smoothness of the value function at c 2 , simply because the first derivative term is integrated 
against the Lebesgue measure which is a diffuse measure. We could also directly use the first part of 
Theorem 6.2 of [16]. □ 

We conclude this section by showing the uniqueness of (c*, c 2 ); recall that the existence was proved 
in Proposition 4.1. 

Proposition 5.3. The maximizer (c^, c 2 ) is unique. 

Proof. Suppose (c*, c 2 ) and (c*, c 2 ) both maximize v CljC2 — c\. We shall show that they must be equal. 
By Lemma 4.5, both [c\, c 2 ) and (cj, c* 2 ) satisfy Case 1 or Case 2 and by Theorem 5.1 we have 

(5.8) v c * )C *(x) = v e * t q(x) = supw 7r (x) x > 0. 

Ti-gn 
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We first show that c 2 = c 2 . Indeed, by Lemma 5.3, v' c * c *(x) < 1 on (c*, c 2 ) and v' & * & *(x) < 1 on (c*, c 2 ) 
while v' c , c *(x) = 1 on (c 2 , oo) and v' e *^*(x) = 1 on (c 2 , oo). Hence if c 2 ^ c 2 , it would contradict with 
(5.8) for the points between c 2 and c 2 . 

In order to show c* = c\, we appeal to the identity f CliC2 (c 2 ) — f C i, C2 ( c i) = c 2 — Ci — /3, which holds 
for any < c\ < c 2 for which w Cl C2 is continuous at c 2 . This together with (5.8) and c 2 = c 2 shows 
^ c *,c*(c*) — f c *,c*(ci) = c* — c*. If 7^ c*, by the mean value theorem, there exists a point between 
these at which v'* „* is one; however, this contradicts Lemma 5.3. This completes the proof. □ 

6. Numerical Examples 

In this section, we confirm the results numerically using the spectrally positive Levy process with i.i.d. 
phase-type distributed jumps [1] of the form 

X t - X = -M + aB t + Y^ Z n , < t < oo, 

n=l 

for some Del and a > 0. Here B = {B t ; t > 0} is a standard Brownian motion, N = {N t ; t > 0} 
is a Poisson process with arrival rate A, and Z = {Z n ; n — 1, 2, . . .} is an i.i.d. sequence of phase- 
type-distributed random variables with representation (m, ot, T); see [1]. These processes are assumed 
mutually independent. Its Laplace exponent (2.1) is then 

^(a) = ds + -a 2 s 2 + A (a(sJ - T) _1 t - l) , 

which is analytic for every s G C except at the eigenvalues of T. Suppose {— £ i>q ; i E X q } is the set of 
the roots of the equality ijj(s) = q with negative real parts, and if these are assumed distinct, then the 
scale function can be written 

(x) = c i,q [ eHq)x ~ e ~^' qX ] and Wiq) ( x ) = Yl Ci 'i [ e * iq)X ~ e_?! ' 9 1 + ~e Hq)x , 

for the case a > and a = 0, respectively for some {C i>q ; i G l q }; see [12]. Here {£ ii9 ; i G X q ] and 
{C7 ij? ; i G X q } are possibly complex-valued. 

In our example, we shall choose a phase-type distribution which does not have a completely monotone 
density. Recall that, in the spectrally negative counterpart [15], the barrier strategy may fail to be optimal 
if the Levy density is not completely monotone. On the other hand, in the dual model, there is no 
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restriction on the Levy measure. We assume m = 6 and 





-5.6546 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 




0.0000 




0.6066 


-5.6847 


0.0000 


0.0166 


0.0089 


5.0526 




0.0007 


T = 


0.2156 
5.6247 


4.3616 
0.0000 


-5.6485 
0.0000 


0.9162 
-5.6786 


0.1424 
0.0000 


0.0126 
0.0000 


, OL = 


0.9961 
0.0000 




0.0107 


0.0000 


0.0000 


5.7247 


-5.7420 


0.0000 




0.0001 




0.0136 


0.0000 


0.0000 


0.0024 


5.7022 


-5.7183 




0.0031 



which give an approximation of the Weibull distribution with density function f(x) = a^ a x a ~ l exp { — {^x) a } 
for a = 2 and 7 = 1, obtained using the EM-algorithm; see [12] regarding the approximation perfor- 
mance of the corresponding scale function. Throughout this section, we let q = 0.05 and let other 
parameters vary so as to see their impacts on the optimal strategy and the value function. 
In our first experiment, we let d = 2, a = or a = 1 with 

Case 1: = 4 and A = 3 
Case 2: (3 = 4 and A = 1 

and obtain the optimal strategies/value functions and confirm the analytical results obtained in the previ- 
ous sections. We choose these parameters so that c* > for Case 1 and c\ — for Case 2. 

Figures 1 and 2 show the results for a = 1 and a = 0, respectively. In both figures, we plot in the 
left column v ClyC2 — C\ with respect to C\ and c 2 and in the right column the value function v c * iC *(-) as a 
function of the initial value x. Recall that the values (c*, 0%) are those that maximize v CliC2 — ci. As can be 
suggested from the contour map of v CliC2 —c±, there exists a unique global maximum and hence Newton's 
method is a reasonable choice of computing the maximizer (c*, d£). For the plots of the value functions, 
the circles indicate the points (c*, v c * t c *(cl)) and (ri^, Vc*^*^)) and the dotted lines the 45-degree lines 
passing through these points. 

In view of these figures, the continuity/smoothness at c\ is readily confirmed; it appears to be differ- 
entiable for the case a = 1 (in other words, the value function is tangent to the 45-degree line) while it 
is continuous for the case a = 0. The non-differentiability for a = is apparent in view of Case 2 in 
Figure 2. At c*, the value function is indeed tangent to the 45-degree line if c\ > 0, while for the case 
c\ = 0, we see that the slope is less than one. These results are consistent with Proposition 4.1. It is 
also confirmed that the slope of v c * >c * is smaller than 1 only at those points inside [tij, cJj], which verifies 
Lemma 5.3 and Proposition 5.2. 

In our second experiment, we take (3 \, and see if the value function converges to the one under 
no-transaction costs as in [8]: 



(6.1) 



-R [q) (a* 



x) 



X, 



if /x > 0, 
if \i < 0, 
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Case 1: (3 = 4 and A = 3 




Case 2: (3 = 4 and A = 1 



Figure 1 . For the case a — 1: (left) v cl>C2 — c\ with respect to c x and c 2 , (right) the value 
function v r * r * as a function of x. 

for any x > 0, with the optimal barrier level 

a . : J(^T 1 (?)>° if ">°- 

[o if M < 0. 

We let A = 3 and consider the case > (by choosing = 2) and also the case /i < (by choosing 
o = 3). 

Figure 3 plots for each case the value function w c * jC *(-) for (3 = 10, 5, 1,0.5, 0.1 (dotted) together 
with the no-transaction case v a *( ) (solid) as in (6.1). The circles on the plots indicate the points 
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Case 1: (3 = 4 and A = 3 
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Case 2: (3 = 4 and A = 1 

Figure 2. For the case with a = 0: (left) v cl:C2 — c\ with respect to ci and c 2 , (right) the 
value function v c * c * as a function of x. 

c 1 ,c 2 

(ci, w c j,c*(ci)), {c*2,v c * ifi * 2 {c* 2 )) and also (a*, v a * (a*)). It is easy to see that the value function is mono- 
tone in (3 (uniformly in x), and converge to the no-transaction cost case as (3 I 0. The convergences of 
both c\ and c* 2 to a* are also observed. In fact, one can prove the convergence of value functions using 
the stability of viscosity solutions. 

Proposition 6.1. Let v 13 denote the value function corresponding to the dividend payment problem when 
the fixed transaction cost is /3 (defined as above), and v the value function when there are no-transaction 
costs. Then converges to v uniformly as (3 I 0. 
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Proof. From the definition of the problem, v" < v and v" is decreasing in (3 and hence it has a point- wise 
limit, which we will call v. The proof is completed if we can show that v is a viscosity super-solution 
of the variational inequality that corresponds to the problem without transaction costs. But this is an 
immediate consequence of the stability result of the viscosity solutions (see e.g. Theorem 6.8 of [19] and 
Theorem 1 of [6]), since we can obtain the variational inequality in the no-transaction case by taking a 
limit in the case with transaction costs. 

To get to uniform convergence from point-wise convergence we just proved, we appeal to Dini's 
theorem to first show it on compacts. This indeed holds because we already know that (v /3 ) and v are 
continuous functions and v 13 j" v as (3 1 0. Now, because the slopes of (v@) and v are all one above c* 2 and 
a*, respectively, and because c\ can be shown to be bounded for any small (3 (thanks to the convergence 
c\ to a* as (3 \, or modifying the proof of Lemma 4.3), the uniform convergence holds. □ 

We also observe in the figures that for fi < 0, c* = 0. This can be shown analytically for any f3 > 0. 
Corollary 6.1. If n < 0, we must have c* = 0. 

Proof. By the nature of the problem the value function f c *,e* is dominated by that of the no-transaction 
cost case. By (6.1), we must have v c * :C *(x) < x for any x > 0. Moreover, because w c *, c *(0) = 0, 
v' c * c * (0+) < 1 and hence in view of Proposition 4.1, we must have c\ = 0. □ 
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